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1 Introduction 



In recent years there has been a growing interest in the calculation of trans- 
port coefficients in quantum field theories at finite temperature. With the 
advent of relativistic heavy-ion colliders, such as RHIC, a proper knowledge 
of transport coefficients has become relevant since hydrodynamical descrip- 
tions of heavy-ion collisions provide a useful tool to analyse the experimental 
data |]l| . In practice the extension of ideal relativistic hydrodynamics to in- 
clude finite transport coefficients is far from straightforward and applications 
to heavy-ion physics have only just begun [|]. 

If the temperature is sufficiently high and the theory is weakly coupled, 
transport coefficients can be computed in a perturbative expansion, employ- 
ing either kinetic theory or field theory using Kubo formulas. It turns out 
that the latter approach requires the summation of an infinite series of higher- 
order diagrams, known as ladder diagrams, which has been a serious draw- 
back for its use. For a scalar theory, the higher-order contributions in the 
loop expansion have been identified and summed in Ref. [Q], using an intri- 
cate diagrammatic analysis, and the leading-order results for the shear and 
bulk viscosities have been found. In Ref. the equivalence of an effective 
Boltzmann equation and the field-theoretical calculation is shown. Using 
a more transparent analysis, the diagrammatic conclusions of Ref. [Q] have 
been confirmed recently [^. So far, the only other transport coefficient for 
which the ladder series has been summed explicitly is the color conductivity 
in QCD P]. The viscosities in the scalar theory and color conductivity in 
QCD have the property that the one-loop contribution and the ladder con- 
tributions are of the same order in the coupling constant. However, for other 
transport coefficients in gauge theories, such as the shear viscosity or the 
electrical conductivity, the ladder contributions are in fact larger than the 
one-loop one Only recently and using kinetic theory, a complete com- 
putation of the leading logarithmic order of these transport coefficients has 
appeared Unfortunately, a full leading-order computation is still lacking. 
Ref. [§] provides a useful guide to the literature. 

Euclidean lattice simulations offer in principle the possibility to compute 
transport coefficients completely nonperturbatively Qj. However, transport 
coefficients are determined by the small frequency limit of zero-momentum 
spectral functions of appropriate composite operators (such as components 
of the energy-momentum tensor) and spectral functions or other real-time 
correlators are not readily available on a euclidean lattice, although recent 
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progress has been made with the Maximal Entropy Method (MEM) pUj, pT[].p| 
For that reason it was proposed in Ref. to introduce instead an ansatz for 
the spectral function and fit the result to the numerical data for the euclidean 
correlator, employing a standard dispersion relation between these two. This 
approach was pursued more recently in Refs. [0, [l^ . 



Motivated by these studies, our goal in this paper is to compute the spec- 
tral function relevant for the shear viscosity at high temperature in weakly- 
coupled scalar (Sec. |^) and nonabelian gauge theories (Sec. H). In the Con- 
clusions we compare our findings with the analysis carried out so far in Refs. 



[P, |I4[ . It is found that the ansatz used in these papers is inadequate and 
we suggest a better one. We also point out a potential problem in the calcu- 
lation of spectral functions at very low frequencies {u 0) from euclidean 
lattice correlators using the MEM approach. 



2 Correlation functions 

We start with a summary of basic relations between transport coefficients, 
spectral functions and euclidean correlators, using the shear viscosity as an 
example. The relations presented in this section are quite general and valid 
for arbitrary transport coefficients. 

The shear viscosity can be defined from a Kubo relation as 

r^=l-\im- f d'xe'-'{[nki{t,^),7rj,i{0,0)]), (1) 
Z[) i^^O w J 

with TTki the traceless part of the spatial energy-momentum tensor. The 
brackets denote the equilibrium expectation value at temperature T. We 
define spectral functions of hermit ian (composite) operators, such as iiki, as 
the expectation value of the commutator, 

PttA^ -y) = {[T^kiix], TTkliy)]), (2) 
and in momentum-space 

pUu;,p) = j d^'xe'^'-'^-^p^^it,^). (3) 

^Note that for classical field theories at finite temperature spectral functions can be 
computed nonperturbatively by numerical simulations directly in real time . 
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Spectral functions obey basic symmetry relation^ P*tt-k{^) = — P7r7r(a^) = 
p^^(— x) and p*^(ct;,p) = p^jrl"^,?) = — PttttI— p) as well as the positiv- 
ity condition u;p^^(co', p) > 0. The shear viscosity is determined by the slope 
at zero frequency: 

1 d 



u)=0 



where Pt,t,{uj) = Pt^t,{u,0). Since the spectral function is odd, we consider 
from now on positive uj only. 

The euclidean correlator (at zero spatial momentum) is given by 

G'i(r) = J d'x{7CMir,^)7TMiO,0))E (r = zt), (5) 

where the imaginary time r G [0, 1/T] and Gf^{l/T — r) = G^^{t). The 
euclidean correlator and the spectral function are related via an integral 
equation, originating from a dispersion relation, 

Gt{r)= r^KiT,uj)pUuj), (6) 
Jo 

with the kernel 

K{t, u) = e'^^n{uj) + e^'^^ [1 + n{uj)] = e'"^^ + 2n{u;) coshcjr, (7) 
obeying K{t,u!) = —K{t, —uj) = K{l/T — t,uj). The Bose distribution is 

n{uo) = -4^- . (8) 

The low-frequency part of the spectral function contains all information 
on the transport coefficient and its effect on the euclidean correlator can be 
estimated quite easily. When u; -C T, the kernel can be expanded as 



UJ I 



l_^T(l-rT) 



+ 0{ujyT'), (9) 



and all except the first term are suppressed. As a consequence we find that 
the contribution to the euclidean correlator from low frequencies, 

Gt^-ir) = 2T r^-^P^, (10) 
Jo ^ 



^ When most manipulations take place in real-space, it can be convenient to define 
spectral functions such that they are real instead of purely imaginary in real-space. Here 
we stick to the usual convention and spectral functions are real in momentum-space. 
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is independent of r. The frequency cutoff c<ja ^ T is introduced here to 
justify the expansion of the kernel. We conclude that the dominant effect 
of the low-frequency region is a constant r-independent contribution to the 
euclidean correlator. 

3 Scalar field 

We consider a one-component massless scalar field with a quartic A0'^/4! 
interaction.^ The one-particle spectral function is 

p{x -y) = m^)A{y)]) =G>{x~y)- G<{x - y). (11) 

The one-particle Wightman functions, G^{x — y) = {(j){x)(j){y)) and (^"^(x — 
y) = {'P{y)4>{^))^ are related to the one-particle spectral function via the 
KMS condition 

G>{p) = [n{p') + 1] p{p), G<{p) = n{p')p{p). (12) 
The traceless part of the spatial energy-momentum tensor reads 

vr,, = d,^d,<P - ^4.5.09.0. (13) 

The lowest-order skeleton diagram that contributes to the spectral function 
for the shear viscosity in Eq. (|^) follows from simple Wick contraction, 

([0^(x), <p'{y)]) = 2 [G>'ix -y)- G<\x - y)] , (14) 

and we find 

p^^iuj) = 1 1 ^ (k ■ \^fn{k')p{k\ k) [p{k' + o;, k) - p{k' - co, k)] , 

(15) 

where we used the KMS conditions (|12D and —n{—uj) = n{uj) + 1. The four 
k's in the integrand arise from the derivatives in rcki- In the remainder of 
this section we compute the one-loop spectral function (|1^) as a function of 
the external frequency. 

•^We assume that the temperature is sufficiently high such that a possible zero- 
temperature mass uiq <C AT^ plays no role. 
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Figure 1: Typical configuration of poles in the complex fc°-plane for the 
evaluation of the one- loop spectral function Pt,t^{uj), (27k) denotes the 
quasiparticle energy (width). 



The spectral function p^^ depends on the one-particle spectral functions 
p, which contain the quasiparticle structure of the theory at finite temper- 
ature. We describe this in some detail since similar, though more compli- 
cated, considerations appear in gauge theories. For hard momenta |k| ~ T, 
excitations are on-shell with energy |k|. For softer momenta screening ef- 
fects become important and hard thermal loop (HTL) resummation yields a 
temperature-dependent plasmon mass, m? = m'^^{l — 3mth/7rT + . . .) with 
m^jj = AT^/24 Finally, collisions in the plasma result in a finite (but nar- 
row) momentum-dependent width 27k ^ m ^ T, changing the one-particle 
spectral function from an on-shell delta function to a Breit-Wigner spectral 
function (see below) |jl6[. 



We may now discuss the one-loop expression ([T5|) . First we note that 
the integral is dominated by hard ~ T momenta. Therefore, for external 
frequencies that are not too small a simple on-shell delta function for the 
one-particle spectral functions suffices to find the dominant contribution. We 
will refer to this region as the high-frequency region. For smaller frequencies, 
however, the arguments of the delta functions come close, producing a so- 
called pinch singularity [0]. The pinch singularity is screened by a finite 
external frequency or width, whichever one is the largest. For very small 
external frequencies, the inclusion of the width^ is essential . This situation 
is sketched in Fig. |l|. The effect of these nearly pinching poles is to enhance 
the spectral function compared to naive estimates. We will refer to the region 
where pinching poles are important as the low-frequency domain. 

We start with the region where the frequency is much larger than the 
thermal width and no pinch-singularity problems are encountered. The one- 



*And of ladder diagrams, see below. 
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Figure 2: Contribution to the spectral function Pt,t,{ijj) / oj^ (full line) from 
decay /creation processes [see Eq. (|17D ] as a function of tu/T, with m/T = 0.1. 
The contribution from the nearly pinching poles in the low-frequency region 
is discussed below. The dashed line shows the contribution proportional to 
the Bose distribution only. The dotted line indicates the asymptotic value. 



particle spectral function can be taken on-shell and is 

Po(A;°,k) = 27re(A;°)(5(A;2 



(16) 



where cuk = vkM-?T? and e{x) the sign-function. HTL effects are included 
in the mass parameter. Evaluating the integrals in Eq. (|15D with the use of 
the delta-functions (the angular integrals are trivial) results in (recall that 
we take u > 0) 



p^^[uj) = 9{uj — 2m) 



ASttcu 



n{uj/2) + 



(17) 



which is shown in Fig. |^. The physical processes are the decay of a zero- 
momentum excitation with energy uj into two on-shell particles with equal 
and opposite momentum, and the inverse process of creation. The decay 
process contributes also at zero temperature and makes the spectral func- 
tion increase as at large frequencies. The threshold at = 2m arises 
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from (simple) HTL resummation. This concludes the analysis of the spectral 
function in the high-frequency domain. 

We continue with the low-frequency region where pinch singularities lead 
to a nontrivial enhancement of the spectral function. We replace the on-shell 
one-particle spectral functions with Breit-Wigner spectral functions: 



27k 27k 



(18) 



The width 27k is determined by the imaginary part of the retarded self energy 

ImSj;(a>k,k) 
= 2^^ ' 

and the dominant contribution at weak coupling arises from two-to-two scat- 



tering from the two- loop setting-sun diagram & Il6[. A convenient way to 



write this damping rate is as |T7| 



7k = 7— 5(|k|/T;m/T), (20) 
c^k 

where 

^- 1536^' ^21) 
determines the parametrical behaviour. The function B contains the non- 



trivial momentum dependence and is related to A defined in Ref. as 
i?(|k|/T; m/T) = (6/7r^)y4(|k|/T; m/T). In the limit of hard momenta and 
small mass [|1^ 

lim lim 5(|k|/T; m/T) = 1. (22) 

For analytical estimates we will neglect the momentum dependence and take 
B = 1. In the results obtained by numerical integration the full momentum 
dependence is included. 

We insert the Breit-Wigner functions into expression (|T5p for the spectral 
function and perform the kP integral by integrating around the poles in the 
complex plane. We preserve only the dominant contributions and discard 
all terms suppressed by (powers of) the coupling constant with respect to 
the leading order contribution. Breit-Wigner spectral functions have four 
poles at complex energy-arguments = ±(0;^ ± i'~fk). From the residue 
of these poles we keep n{uji^ ± Z7k) ~ n{uJk). The Bose distribution n{k^) 
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has poles along the imaginary axis at = 27rinT, n G Z. However, the 
residues at these poles are subdominant compared to those from the poles 
of the Breit-Wigner functions. Hence we do not include these contributions. 
After performing the integral we find 



with 



and 



d^k 



(2vr) 



3' 



(23) 
(24) 



A^>k) = - (25) 

For sufficiently large u and in the limit of small coupling (width) we may 
take /(tu, k) 47r(5(c<j — 2a;k) and the result of the previous calculation in the 
high-frequency region is recovered, as it should. 

Let's now consider the low- frequency region, < m. Here we may ap- 
proximate /(ct;,k) ^ 27k/co'k and expand the difference between the Bose 
distributions to find 

pU^) = -| / K^'M 2 (0 < ^ < ^)- (26) 

We note that the last factor controls the pinch singularity in precise agree- 
ment with Fig. |I[ Although in principle p,r7r('^)/r^ may depend on the three 
dimensionless combinations uj/T, 7/T, and m/T, we find that in practice 
it only depends on u/'-f and m/T. Since the integral is dominated by hard 
momenta the m/T dependence is subdominant and the natural parameter on 
which the spectral function depends in this region is 1^7/7. It is easy to see that 
the spectral function has a local maximum at u; ~ 7 and p-rrni^^ ~ ~ 1- 

For very small frequencies ^ 7 we expand and obtain 

p^Au^)/T' = ar(-)+^(-X + ... (0 < c « 7), (27) 

with 



7/ 3! \7 



2 /• |k|S, ,7 5!C(5) 
ai = --— / —n M— ^ 28 
3r4 Jk 7k Svr^ 

,.J,/^„V.)m^-^. (29) 



T4 ul \^^J 27T 
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Figure 3: Ladder diagrams that contribute to p-Kiriyj) in the scalar theory. In 
the low-frequency region, < 7, the pinching-pole contributions from the 
ladder diagrams are equally important as the one-loop contribution. When 
^ 7, pinching-pole contributions from ladder diagrams are suppressed. 



where the ~ indicates that the final integrals are evaluated by neglecting the 
remaining momentum dependence of the damping rate (i.e. taking i? = 1) as 
well as the thermal mass (the first approximation has numerically the largest 
effect). From this the one- loop viscosity follows as 

1 f Ikl^ , 1 2C(5)T^ , , 

30 Jk 7k vr^ 7 

However, as is well-known these one-loop results are not complete and 
a ladder summation is required to obtain the complete leading-order result 
(see Fig. ^). Due to nearly pinching poles, each additional rung in the ladder 
contributes with a factor A^T/7 ~ 1 and is therefore not suppressed. The 
effect of ladder summation is to change the coefficients ai, 03, . . ., but not the 
parametric dependence on the coupling constant. 
In the region 7 ^ < m we find 

/ 8 /■ |k|^ ,1 ^7k 8C(3)7 I < \ 

PMuj)/T =-—^ — r2(^k)— ^ — 2 (7<^<m). (31) 

01 Jk "-"^k TT U! 

In this frequency interval pinching-pole contributions from ladders are sub- 
dominant since each additional rung comes with a factor X^T/uo ^ 1. The 
perturbative part of the three-loop ladder diagram (i.e. with a single rung) 
also contributes at this order and has the same }?T jio behaviour as we find 
above.0 Fig. ^ shows the contribution to the spectral function from the nearly 
pinching poles in the low-frequency interval, obtained by numerical integra- 
tion of Eq. (|26D with the full momentum and mass dependence. Note that the 
natural dimensionless combinations in the low-frequency region differ from 
those in the high-frequency region (compare Figs. |^ and 

^We thank Guy Moore for pointing this out. 
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Figure 4: Contribution to the spectral function Pt^t^{uj)/T^ (full line) from the 
nearly pinching poles in the low- frequency region as a function of uj/'-f, ob- 
tained by numerical integration of Eq. (|26|). The dashed line is the analytical 
result (|3T|) when 7 -C < m, neglecting nontrivial momentum dependence 
and finite mass corrections. The inset shows a blowup. The viscosity is 
determined by the slope for a; 0. 

We may now combine the results obtained so far to construct the complete 
one- loop spectral function at high temperature in the weak-coupling limit. 
The spectral function can be written as the sum of the contributions discussed 
above: 

P..M=plTM + p'fM, (32) 

where p^°^ represents the contribution from the nearly pinching poles in Eq. 
([2B|), dominating at low frequencies, and pjj'^^ is the contribution from de- 
cay/creation processes in Eq. (0), dominating at higher frequencies. 

We find that for very small frequencies the spectral function rises quickly 
as w/'j. It reaches a maximum of order 1 (in units of temperature) at ~ 7 
and decays then slowly as 7/0;. Around the thermal mass the contribution 
from decay/ creation processes enter and the spectral function increases again. 
Note that the results obtained in both frequency domains smoothly match 
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Figure 5: Complete one-loop spectral function Pt,t,{ijj) / uj^ (full line) as a 
function of cu/T, with m/T = 0.1. The dashed line is the contribution from 
decay/creation processes and vanishes below uj = 2m, the dot-dashed line is 
the contribution due to the nearly pinching poles at lower frequencies. The 
dotted line indicates the asymptotic value. 



parametrically at a; ~ m, since Pmri^^ ~ 3m) /T^ ~ A^/^ both from the 
low- and the high-frequency calculation. For large u the spectral function 
increases as u;^, due to the zero-temperature decay process. Ladder diagrams 
do not affect this characteristic shape. In Fig. ^ we present the complete one- 
loop spectral function as a function of cu/T.^ In order to combine the low- 
and the high-frequency contribution in one figure, we show Pt^t^{uj) / uo^ , which 
enhances the contribution at lower frequencies. 

We are now ready to calculate the euclidean correlator using Eq. (H). 
Because G^^ij) depends linearly on the spectral function, we write it as a 
sum of two contributions, G^^ = G^^°^ + Gf^^^^, and discuss each term 
separately. We start with the contribution due to decay /creation processes 

We used here that m/T = 0.1 corresponds to A ~ 0.267 and j/T ~ 1.48 • 10"^. 
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which reads 



K{t,uj)- 



2m 



2ti 



n{u/2) + 



(33) 



It is easy to see that the mass plays only a sub dominant role and finite-mass 
corrections are suppressed by m^/T^. Therefore we take m = which yields 



G 



(r) 



2^5 



967r2 



n[x) 



n{x/2) + 



(34) 



where x = uj/T and s = rT. The remaining integral can be performed and 
we find 



3 sin^ u 



T 



{(tt — u) [11 cosM + cos3m] + 6sinu + 2sin3M} , (35) 



where u = 2'ktT . An approximate but illuminating expression for the eu- 
clidean correlator can be obtained by noticing that the integral in Eq. (|3^ 
is dominated by hard frequencies, x > 1/s with < s < 1, such that the 
Bose distributions may be approximated with Maxwellian ones, n{x) ~ e~^. 
In that case we find 



87r2 



■ 1 

7^ + 



1 



{l/T-rf (3/2T- 



r 



(l/2T + r) 



• (36) 



This approximate expression differs less than 2% from the exact result (pSj). 
The dominant 1/r^ behaviour of the correlator arises from decay at zero 
temperature. Finite temperature is manifested mainly through the reflection 
symmetry, r — >■ 1/T — r. At the central point tT = 1/2, we flnd the 
contribution to the euclidean correlator from the decay/creation processes to 
be 



47r^ 

Is 



1 - 



25 

8^r2 



+ 



(37) 



The effect of a flnite mass is to lower the mimimal value of the correlator. 

The contribution to the euclidean correlator from the nearly pinching 
poles in the low-frequency region can be found easily from Eqs. (^ and (p6[ ) 
by interchanging frequency and momentum integrals. Since p^°^ gives the 
dominant contribution to the spectral function up to frequencies of order m. 
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Figure 6: Logarithm of the euchdean correlator G^^{t) /T^ as a function of 
tT, obtained by numerical integration of Eq. @ with (p2D, for m/T = 0.1 
(full line). The dashed line represents the analytical expression for m = 0, Eq. 
([39|); it cannot be distinguished from the full result, except in the inset. The 
dot-dashed line shows the contribution from decay/creation processes only, 
obtained by numerical integration of Eq. (|33|) . The inset shows a blowup 
around tT = 0.5. 



whereas p^'^^^ dominates for higher frequencies, we can use the expansion 
for the kernel and obtain 



G 



E,low , 



r) 



k ^k 



— ri(^k) / — 







27r uj uj^ + 4% 



SvrJk ^k V27ky 



47r2 

Is 



rp5 



25 m2 
8^T2 



+ . . . 



(38) 



with c<ja ~ m. The error that is introduced by expanding the kernel is negli- 
gible (see Fig. ^). We note that in our one- loop calculation this result is at 
leading order independent of the coupling constant. Similarly, while ladder 
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diagrams determine the precise shape of the spectral function, the effect on 
the euchdean correlator G^^'°™(r) appears only in subleading corrections in 
the coupling constant, as can be seen with the kinetic approach []18[. There- 



fore the low-frequency contribution to the euclidean correlator is constant 
and of order one (in the appropriate units) and insensitive to details of the 
ladder summation. 

Combining the low- and high-frequency contribution to the euclidean cor- 
relator at high temperature and weak coupling we find, to leading order in 
the coupling constant, 

G^At) = E — {(tt — u) [11 cosu + COS 3m] + 6 sinw + 2 sin 3u} + — — T^, 

3 sin u 45 

(39) 

with u = 2t{tT . Corrections due to the finite thermal mass are suppressed 
by m? /T"^. The euclidean correlator is minimal at tT = 1/2, and G^^(1/2T) 
receives contributions of the same order from both the high- and the low- 
frequency region (at leading order in the coupling they are equal). A com- 
parison between the analytical result at m = and the full result obtained 
by numerical integration of Eq. (P) with Eq. (|32D in the presence of a finite 
mass is shown in Fig. ||. 



4 Nonabelian gauge fields 

We leave the scalar case and consider nonabelian SU{Nc) gauge theory. The 
traceless spatial part of the energy-momentum tensor is 



vr., = FfF^^ - -S,,F^'^F^^. (40) 



1 

— ( 
3 

The coupling vertex between the operator TCij and two gluons with incoming 
momenta P,K and indices (/x, a), (z/, 6) respectively can be read from (pO]) 
and we find^ 



5^iu [ Pikj + Pjki - ^(5ijp ■ k 



+P ■ K(^6ifj,6ju + 6iuSjfj, - -SijSk^,Ski,^ - {piK^5ju + PjK^6i, 



^ In order to arrive at the basic one-loop expression (Q) below we use here the 
imaginary-time formalism with P = (p4,p), the Matsubara frequency w„ = — p4 — 2'KnT 
[n £ Z) and P ■ K = p^k^ + p • k. 
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Figure 7: Diagrams that contribute to the spectral function Ptttt at higher 
order. These diagrams are special for a nonabelian theory. 



(41) 



In the nonabelian theory rcij also couples to three and four gluons, which 
leads to diagrams as depicted in Fig. |^. However, these contributions are 
suppressed by powers of the coupling constant and will not be considered 
further. 

We use the Coloumb gauge in which the gluon propagator reads 

D,uiP) = PJ,(p)At(P) + 6,,6,,ALiP), (42) 

with Vfj{p) = Sij — pipj, Vj^ = PJ4 = 0. The transverse and longitudinal 
components have the following spectral representations 

A,(p)=- r^^ii;^^. A,(P)=±+ r^^^ii^. (43) 

2% lUJn -UJ 2-K lUJn - UJ 

The one-loop contribution to the spectral function reads, after evaluating the 
Matsubara sum and taking all frequencies real again, 

p..(^) = ^ I 0-^[n{k')-n{k° + u;)]S^V,{k,u)pT{k',k)pT{k° + u;,k) 
+V^2(A;)pt(A;°, k)pi(A;° + u,k) + V3{k)pL{k\ k)pi(A;° + u, k) 1,(44) 



where = — 1 is the number of gluons and Vi{k, uj) = 7|k|^ — 10k^A;''(A;'' + 
cj) + 7(fc0)2(fc0 + cj)2, V2ik) = 6k2(fc0)2 and Vs{k) = -32|k|^ This expression 
is the equivalent of Eq. (^) in the scalar theory. 

As in the scalar case, we start the analysis for external frequencies lj 
that are sufficiently large such that no pinch singularities are present. The 
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collective (HTL) effects in a nonabelian plasma can be incorporated in 



the one-particle spectral functions, which are however more complex than in 
the scalar case. The spectral function for transverse gluons can be written 

as m.m 



PT{k\ k) = 27rZT(|k|) {5[k' - uorm)] - + MM)]] + /3t(A:°, |k|). 

(45) 

The delta functions describe propagating quasiparticles with a dispersion 
relation ut and a residue 

7 = Mk)[ujUk) - k'] .... 
' 3ujl,ujUk) - [u^k] - k^Y ' ^ > 

where = is the plasma frequency squared. For small and large 

spatial momentum one finds 

6 
5 

uj^{k) ~ P + m^, Zrik) ~ 1/(2A;) {k ^ oo), (48) 

where = |a;pi is the asymptotic gluon mass squared. The (3t function 
describes Landau damping and is nonzero below the light-cone only. The 
spectral function of longitudinal gluons Pl(^°, k) has a similar form |2^, 21 



ujUk)^coli + -e, Zt(A;) ~ l/(2cupi) (A; 0), (47) 



We start studying the contribution to Eq. (0) when both gluons are 
transverse. There are three parts, depending on whether we take from the 
one-particle spectral functions the delta functions of the quasiparticles (pole- 
contribution) or the P function of the Landau damping (cut-contribution). 
When both gluons are quasiparticles the integrals can be done with the help 
of the delta functions and we find the pole-pole contribution for transverse 
gluons 



(49) 

Here we use the notation /o = f{uj/2), with f{u) defined as the inverse of 
the transverse dispersion relation, f[uT{k)] = k. The function f{u) vanishes 
when u < tUpi, and 



/(«)^^^(«2-^2^) (u^oo,,), (50) 

3^pi 

/(n)~M-^ (u^oo). (51) 



17 



Again, as for the scalar case, the HTL resummation produces the threshold 
in the spectral function (^9[) for soft external frequency. We find that 

pZ{u^ - gT) - g'T\ (52) 

For large uj the spectral function behaves as 

I 



n(o;/2) , 2 



(u;>a;pi), (53) 



which is, up to the prefactor, what we found in the scalar well. 

The contribution when one transverse gluon is a quasiparticle and the 
other undergoes Landau damping (pole-cut contribution) is 

= ^ Tdu -4^,ZT[f{u)Mu -u)- n{u)mu - UJ, f{u)] 
X [7f\u) - lOuiu - uj)f\u) + 7m2(m - ujf] . (54) 

For soft UJ ~ gT the dominant contribution arises when the energy u is hard, 
and we may use 

(3^[u~ujJ{u)]^^-^% (55) 
as well as other simplifications given above. We find 

plliuj ^ gT) -^ujI, / du u'n'iu) ~ g'T\ (56) 

For soft external frequencies the pole-cut contribution dominates over the 
pole-pole contribution (|5^. 

For hard frequencies ~ T, the (3t function determines the lower inte- 
gration limit to be uo = uj/2 + 3a;pi/(4ci;). As a result, u is always hard and 
we can simplify the integrand to arrive at 



n{u — uj) — n{u)]/3T{u — uJ, u). 



pll{uj ^T) = ^ duu' {Au' - Auuj + 7oj') 

(57) 

It is convenient to substitute u = uj{z + 1)/2 such that 

(58) 
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where z± = {z ± l)/2 and a = SUp^/ {2uj'^) . The dominant contribution 
comes from the lower integration hmit {z a) where the /^y function can 
be approximated as 

Pt{ujz_,ujz+) - (59) 

Using this expression we find that 

plUoo-^T)^g'THn{l/g). (60) 

For hard frequencies the pole-cut contribution is therefore suppressed com- 
pared to the pole-pole contribution (0). Finally, we found that the remain- 
ing cut-cut contribution when both gluons are transverse is suppressed with 
respect to the pole-pole contribution when uj is hard and to the pole-cut 
contribution when u is soft. 

In a similar way we have analysed the remaining longitudinal-transverse 
and longitudinal-longitudinal contributions in Eq. (0) with the result that 
they do not modify the conclusions drawn from the transverse-transverse 
contribution analysed above. In particular, for hard frequencies u ^ uipi 
transverse gluons dominate and the spectral function is given by Eq. ( p3D . We 
do not need to be more explicit about the region where u < gT because the 
dominant contribution in this region arises from the pinching poles, screened 
by a finite width and/or external frequency, as we will show now. 

For small external frequencies < u; < gT the loop integral is dominated 
by the region of hard momentum and we only need to consider transverse glu- 
ons, since the residue for longitudinal gluons vanishes exponentially. In order 
to avoid pinch singularities we follow the same steps as in the scalar case and 
substitute for the one-particle spectral function a Breit-Wigner function, see 
Eq. (0). For transverse on-shell gluons with hard momentum the disper- 
sion relation is = ^/k^ + and the leading (momentum-independent) 
contribution to the damping rate is |^| 



j = ^N,T\n{l/g), (61) 

where the logarithm is sensitive to the magnetic mass, ln(ci;pi/minag) ~ lii(l/5f) 
with minag ~ The mass rrioo only plays a subdominant role and is ne- 

^In QED there is no magnetic mass which could regularize the logarithmic divergence 
of the leading contribution to the damping rate. It turns out that the electron damping 
rate is ill defined lESl. 
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glected below. Evaluating the integral over exactly as in the scalar case, 
we find 



= -y/ ^{ Klk|) -n(|k| -c.)]/(c.,k)A(c.,k)-l^L_^ 

-[u^-uj]], (62) 



where I{ui, k) was defined in Eq. ( p5|) and A{uj, k) = k^ (4k^ — 4c<j|k| + 7uj'^). 
For < Upi this expression simplifies to 

= / kV(|k|)^^ (0 < < ..pi), (63) 

which is 2(i^ times the scalar result. In the region 7 ^ < Upi we find 

p^^^ /T =-— ^- / kn(k) = — — d^-. 64 

As in the scalar case, the three- loop ladder diagram with one rung contributes 
in this region at the same order. Note that the contribution from the pinching 
poles at u; ~ gT, 

pU^^gT)r^gTHn{l/g) (65) 

actually dominates over the contribution ( ^6]) from the collective (HTL) ex- 
citations in this region. The shear viscosity in the one-loop approximation 
follows from (j6^) as 

Vl-loo, = (66) 

However, for the complete spectral function at small frequencies < < 7 
the effects of ladder diagrams must be taken into account. As was men- 
tioned in the Introduction, kinetic theory predicts that the shear viscosity 
is parametrically larger than the one-loop result ( pBD [[^, ^ There- 
fore, for vanishing frequency the spectral function due to ladder diagrams 
is expected to be larger than the one-loop contribution and to behave as 

gladder (^^ 0) = 20// with § 

- 1 
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Hence the slope of the spectral function is much steeper close to the origin, 
compared to the one-loop result. Up to frequencies ci; ~ 7 the actual be- 
haviour of the spectral function depends on the pinching-poles contribution 
of ladders diagrams. When 7 <C u; < ujpi pinching-poles contributions of 
the ladders are subdominant by a factor g'^T/uj and the spectral function 
decreases as g^T jio until it meets with the rising contribution from Eq. (^31). 
For large lo the spectral function increases as cu^. 

The euclidean correlator at high temperature can be computed as in the 
scalar case and we find, at leading order in 

G [Tj/T = E — [in — u) 11 cosM + cos 3m + 6 sinw + 2 sin 3m|H , 

sin u 45 

(68) 

with u = 2'ktT . The last constant term reflects the low-frequency region. 



5 Conclusions 

We have studied the spectral function relevant for the shear viscosity in scalar 
and nonabelian gauge theories at high temperature as a function of the ex- 
ternal frequency. While for small frequencies ladder diagrams are important 
in the scalar case and essential in the nonabelian case, for higher frequencies 
a simple one-loop computation yields the dominant contribution. 

We found that the spectral function has a characteristic shape: for small 
frequencies the spectral function rises, reaches a local maximum and de- 
creases as l/oj. This contribution is due to scattering processes in the plasma 
and is enhanced due to nearly pinching poles. We referred to this contribution 
as the low- frequency contribution. For higher frequencies, decay/creation 
processes dominate and the spectral function increases essentially as uj^. This 
contribution is referred to as the high-frequency contribution. 

In order to extract transport coefficients from euclidean lattice correla- 
tors, a simple ansatz, essentially a Breit-Wigner spectral function, was intro- 
duced in Ref. to model spectral functions of components of the energy- 
momentum tensor. The resulting three-parameter ansatz was subsequently 
used in Refs. |n| to determine transport coefficients in hot gauge theo- 



ries from lattice simulations. Unfortunately, as we have seen in this paper, at 
high temperature spectral functions of composite operators do not resemble 
simple Breit-Wigner functions at all. In order to improve this analysis, we 
propose therefore to use a different ansatz, which is written as the sum of 
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two terms: 

P..H=ptTM + p'f M- (69) 

The high-frequency part can be described by Eqs. (^) or ([TTD with m as a 
possible free parameter. For the low-frequency part we note that the spectral 
function is odd, increases linearly with uj for small uj and decreases with 1/uj 
for larger uj. A simple ansatz reflecting this is 

p^^[x)T^ = x- , x = ujT, (70) 

with bi = Ci = 0, i > n for given n. The viscosity is given by rj/T"^ = bi/20. 
The ansatz (|69|), with free parameters m, 6j, and q, should be used in Eq. @ 
to fit the corresponding euclidean correlator to the numerical results. When 
insisting on a three-parameter fit, one may take n = 1 which leaves m, bi 
and Ci to be determined. 

Concerning the euclidean correlator, we found that the dominant r de- 
pendence is determined by the high-frequency part {u > T) of the spectral 
function. However, around tT = 1/2 both the high- and the low- frequency 
regions of the spectral function contribute at the same order. The low- 
frequency contribution is of special importance since transport coefficients 
are determined by the slope of the spectral function at zero frequency and 
a precise calculation of the spectral function at low frequencies is therefore 
essential. We found that no matter how complicated the spectral function 
up to frequencies u ~ gT {g -C 1) might be, its contribution to the eu- 
clidean correlator will be a r independent constant. This latter feature poses 
a severe challenge for the MEM approach, since the spectral function at low 
frequencies should be reconstructed from the knowledge of a single constant 
alone. It turns out that this constant carries information on the transport 
coefficient (or more generally on ladder diagrams) only in subleading con- 
tributions. As a result euclidean correlators are remarkably insensitive to 
transport coefficients, which makes it extremely difficult to extract those at 
weak coupling. 

We emphasize that these results are not specific for the correlator we 
considered here. For instance, in the first paper of Ref. ||10[ the current- 
current correlator relevant for thermal dilepton rates in QCD was studied 
on the lattice and an enhancement in the central value of the euclidean 
correlator compared to the free result was observed. This enhancement might 
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be accounted for by the pinching-poles contribution in the low-frequency 
region of the spectral function. ^ 
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